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GENERIC AND NON-GENERIC BEHAVIOR OF SOLUTIONS TO 
DEFOCUSING ENERGY CRITICAL WAVE EQUATION WITH 
POTENTIAL IN THE RADIAL CASE 

HAO JIA, BAOPING LIU, WILHELM SCHLAG, GUIXIANG XU 


Abstract. In this paper, we continue our study m on the long time dynamics of radial 
solutions to defocusing energy critical wave equation with a trapping radial potential in 
3 + 1 dimensions. For generic radial potentials (in the topological sense) there are only 
finitely many steady states which might be either stable or unstable. We first observe that 
there can be stable excited states (i.e., a steady state which is not the ground state) if the 
potential is large and attractive, although all small excited states are unstable. We prove 
that the set of initial data for which solutions scatter to any one unstable excited state 
forms a finite co-dimensional connected manifold in energy space. This amounts to the 
construction of the global path-connected, and unique, center-stable manifold associated 
with, but not necessarily close to, any unstable steady state. In particular, the set of data 
for which solutions scatter to unstable states has empty interior in the energy space, and 
generic radial solutions scatter to one of the stable steady states. Our main tools are (i) 
near any given finite energy radial initial data (uo,ui) for which the solution u(t) scatters 
to some unstable steady state (j) we construct a manifold containing (mo,mi) with the 
property that any solution starting on the manifold scatters to (j); moreover, any solution 
remaining near the manifold for all positive times lies on the manifold (ii) an exterior energy 
inequality from OUHlIIl]. The latter is used to obtain a result in the spirit of the one-pass 
theorem [22], albeit with completely different techniques. 


1. Introduction 


Fix f3 > 2. Define 

Y := I V e C'(M^) : V radial and snp(l + |a;|)^|U(a;)| < cxo [■, 


with the natural norm 

||U||y := sup(l + |x|)^|U(x)|. 

xGK® 


We study solutions to 


dttu — Am — Vu + = 0, 


( 1 . 1 ) 


with initial data 'S’(O) = (mqjWi) G x Since for a short time the term Vu can 

be considered as a small perturbation, by adaptations of results in [HUllllTSlES] we know for 
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any initial data {uq,U]) G x L^(R^), there exists a unique solution 

u{t) e ^([0, oo), H^) n T) X 

for any T < oo to equation fll.ll) . Moreover, the energy 

■|Vn|2 {dtuf Vu^ 


£(V(«)) := 




x, t) dx 


is constant for all time. Our main interest in this work is the long time behavior of u{t) 
under radial symmetry. If the operator —A — V has negative eigenvalues, then equation 
fll.ip admits a nontrivial ground state Q > 0, which is the global minimizer of 




|V0|2 1/(^2 ^ 




... 2 2 

It has negative energy. The linearized operator around Q 

Cq := -A-V + 5Q^ 


dx. 


has no negative or zero eigenvalues, and no zero resonance. Consequently by well-known 
dispersive estimates for Cq we know Q is asymptotically stable: solutions with initial data 
close to {Q, 0) will scatter to {Q, 0). We remark that in our work by generic choice of poten¬ 
tials, all steady states are hyperbolic 0 and consequently spectral stability implies asymptotic 
stability by well-known dispersive estimates for the associated linearized operator. Hence 
we will not distinguish the two notions below. In addition to the ground states Q and —Q, 
there can be a number of “excited states” with higher energies (see Appendix A of ra). 
which are changing sign steady states to equation fll.ip . Surprisingly, some of the excited 
states can be stable as well, although all small excited states can be shown to be unstable 
(see Section [2] below). More precisely we have the following result. 


Theorem 1.1. There exists an open set O C V, such that for V ^ O, there exists an excited 
state 0 to equation M.l\) which is stable. 


Roughly speaking, this is due to the stabilizing effect of the nonlinearity as a result of its 
defocusing nature, and the instability is mainly due to the potential. Hence if the excited 
state is large, the nonlinear stabilizing effect may dominate and the resulting dynamics 
around that excited state could become stable. 

Due to the presence of many steady states, in general the global dynamics can be quite 
complicated, even in the radially symmetric setting. [16] establishes the following result 
characterizing the long time dynamics of radial hnite energy solutions. 

Theorem 1.2. Let (uo,Ui) G x be radial. Denote 

E = {(0, 0)1 (0, 0) is a radial steady state solution to equation fll.ip }. (1.2) 

Let u G C{[0,oo),H^)r\L^Ll^{[0,T) x M^) for any T < oo be the unique solution to equation 
M.l\) with initial data (m( 0), (9tM(0)) = {uq,Ui). Then for some radial finite energy solution 


^This means that the linearized operator around the steady state has neither zero eigenvalues nor zero 


resonance. 
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{u^,dtu^) to the linear wave equation without potential^ 

(LW) duU-Au = 0, 

we have 

lim inf \\{u{t),dtu{t)) - - {u^(t),dtu^{t))\\^i ^^2 = ^. (1.3) 

t^oo ( 0 ,o)es 

Moreover, for V in a dense open set Q, (ZY, there are only finitely many radial steady states 
to equation In this case, there exist a steady state solution (0, 0) and some solution 

{u^,dtu^) to the linear wave equation without potential, such that 

lim \\{u{t),dtu{t)) - (0,0) - {u^{t),dtu^{t))\\pi ^^2 = 0. (1.4) 

t^OO 

We remark that we can in fact choose the set C F such that for any 1/ G 11, all steady 
states are hyperbolic!^. We £x this choice of H below. Theorem ll.2l is a particular instance of 
the soliton resolution conjecture for general dispersive equations, which has been intensively 
studied for many dispersive equations. We refer the reader to [Sl[IU] and references therein 
for results on the focusing energy critical wave equation, [SllGlISlfTSlfT^ and references therein 
for results on equivariant wave maps, and [23[29] for results on Schrodinger equation with 
potential. The difference between these works on this lies with the defocusing nature of our 
equation which precludes any blowup. In other words, the flow on phase space is global in 
time, and together with [16] the present work establishes a complete description of the long 
term dynamics as well as a decomposition of the global data set into components which lead 
to distinct hnal states. 

The result in [T 6 | proves convergence for all radial solutions, thus establishing the so-called 
soliton resolution in the setting of equation fll.ll) . The proof relies crucially on the channel of 
energy inequalities for the linear wave equation, introduced in the works of Duyckaerts, Kenig 
and Merle [9l[in|. This tool implies, amongst other properties, that all non-stationary radial 
solutions emit a positive amount of energy into large distances (the “far held”). The main 
local decay mechanism for equation fll.ll) is the dispersion of energy into large distances, and 
the channel of energy inequalities provide a powerful tool to quantify such effects. In fact, 
due to the presence of the potential which destroys many of the favorable algebraic identities 
of virial type, 0 the channel of energy inequality is perhaps the only tool currently available to 
measure dispersion in this context. As a consequence, in absence of radial symmetry, where 
the channel of energy inequalities (see m) become less effective, we have little knowledge 
of the “compact solutions”, i.e., solution u(t) with the property that {u(t),t G M} is pre¬ 
compact in X L^. This should be compared to the focusing energy critical wave equation, 
for which one knows (albeit along a sequence of times) modulo symmetries, that compact 
solutions converge to some steady state. @ 

^We often call such linear solutions free radiation. 

^In Theorem 6.1 |16j . we only showed that in the radially symmetric case, the linearized operator has 
neither a zero eigenvalue nor a zero resonance when restricted to radial functions. This leaves the possibility 
of having zero eigenvalue or zero resonance when we consider nonradial functions. We will address this issue 
in Section 2. 

^The virial type identities can still be of some use even in this context, see [27]. 

recent result of Duyckaerts, Kenig and Merle [12] . under certain non-degeneracy assumptions, com¬ 
pletely characterizes all compact solutions as Lorentz transformations of steady states. 
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In this paper our main goal is to obtain refined descriptions of the global dynamics of 
solutions to equation (11.11) in the radial case. Let us denote 

X := {(wo,Mi) e X : {uo,ui) radialj . 

We establish the following result. 

Theorem 1.3. Let 12 be an open dense subset of Y such that equation fO) has only 
finitely many steady states, which are all hyperbolic, and let S be the set of radial steady 
states. Denote lt{t) := {t){uo,Ui) as the solution to equation M.l\) with radial initial data 
{uo,Ui) G i/rad ^ -^rad(®^)- (0,0) G S, define 

Mtf, := |(mo,Mi) G x : '^{t){uo,Ui) scatters to (0,0) as t ^ +cxo| . (1.5) 

Denote 

:= -A - V + (1.6) 

as the linearized operator around fi. If has no negative eigenvalues, then is an open 
set C X L^^^(M^). If Cfj , restricted to radial functions has n negative eigenvalues, then 
is a path connected manifold C hf^ad ^ -^rad(®^) codimension n . 

Remark. This result shows that each unstable excited steady state attracts a finite co- 
dimensional manifold of solutions, hence scattering to unstable excited states is non-generic. 
If has no negative eigenvalues, then 0 is stable. This is a relatively straightforward 
consequence of the known dispersive estimates for (see a)0 and standard perturbation 
arguments. 

On the other hand, if has negative eigenvalues then the local dynamics near (0, 0) is 
nontrivial. Thanks to [21] and reference therein, it is now well-known in a small neighborhood 
(in the energy space) of (0, 0) that we can construct a center-stable manifold, on which the 
solution scatters to (0,0). Off that manifold, the solution will exit the small neighborhood 
in finite time. In particular, this center-stable manifold is unique. Generally speaking, 
after exiting the small neighborhood, we lose control on the dynamics based on perturbative 
arguments alone and some global information is needed. While the one-pass theorem 
provided this global information in [21], here it is the channel of energy inequality that 
allows for the key global control on the solution after the exit time. We will provide further 
explanations below. 

Let us briefly outline the main ideas in the proof of Theorem 11.31 Take any unstable 
steady state (0,0) G S and a radial finite energy solution lt{t) with initial data {uq,Ui) 
which scatters to (0,0), i.e., for some radial solution lt^(t) to the linear wave equation 
(LW), we have 

||5/(t) - lt^{t) - (0,O)||^i^^2(R3) = 0. (1.7) 

We first show in a small neighborhood Bfi{uo,Ui)) C II^i,^^ x L^^j(M^) there exists a local 
manifold Ai, such that any solution lt{t) with initial data on this manifold remains close to 

®Due to our relatively mild decay assumption on V, the dispersive estimates we need are close to optimal, 
hence the need for the work [i] which requires less decay on the potential than, say, |L(a;) | < (i_|_ | ;(.|p+ that 
is usually required for the boundedness of wave operators in some other works. 
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it(t) in X for all positive times and also scatters to (0,0). Moreover, this manifold has 
the following uniqueness property: any radial hnite energy solution it {t) which stays close to 
lt{t) for all positive times necessarily emanates from A4. The construction of this manifold 
differs from the usual ones in that this is not a center-stable manifold around a steady state. 
In fact, since the energy of the solution lt{t) may be much higher than that of (0,0), the 
free radiation it^ may contain a large amount of excess energy. One new technical aspect 
is that in addition to using fll.7p . we also need the space-time control on the radiation term, 
such as 

u-(^e ^Ll^{[0,oo) xR^). (1.8) 

fll.Sp is of course expected, but was not usually mentioned in the literature. With the help of 
fll.Sp . the construction of Ai follows from standard techniques. The next step is to describe 
the dynamics of solutions starting in Ui)) C ^ , but off the manifold Ai. 

This is where we need the global control provided by the channel of energy inequalities. Take 
any solution it(t) starting in i?e((Mo,Mi)) C X L'^^^(Rt) and off the manifold (possibly 
with a smaller e), then by the property oi Ai, u{t) — it(t) will have energy of a hxed size 
at some time t, no matter how small ^(0) — ^(0) is. We will show from this that it(t) will 
emit a hxed amount more energy than it, thanks to the channel of energy inequality. The 
main difficulty is that since it (t) may have already emitted a large amount of energy in order 
to settle down to 0, we need to distinguish the new radiation from the old radiation. This 
is done with careful perturbation arguments as follows. Choose (uoWi) very close to {uo,Ui) 
so that the solutions it(t) and it (t) remain close for a sufficiently long time. During this 
time, the radiation has propagated sufficiently far from the origin (with the bulk of energy 
traveling at speed ~ 1). In the hnite region, the solution it is just a small perturbation of 
(0,0). Due to the assumption (uoWi) t after another long time it(t) will deviate from 
(0, 0) in the hnite region by a hxed amount. Then we apply the channel of energy inequality 
to show that it {t) emits a second piece of radiation, which is supported very far away from 
the hrst radiation. Hence, in total it (t) emits quantitatively more energy into spatial inhnity 
although the energy of it (t) can be chosen arbitrarily close to that of lt{t). Consequently, 
it [t) has less energy than (0, 0) in the hnite region for large times, and must scatter instead 
to a steady state of lower energy, not (0, 0). This establishes the proposition that in a small 
neighborhood of (mo,mi), only initial data on Ai can lead to solutions scattering to (0,0). 
Thus Ai^ is truly a global manifold in x L‘^^^(R^) whence Theorem 11.31 The fact that 
Ai^is path connected follows from a perturbation argument which we present at the end of 
Section [5l 

1.1. Some open questions. Our investigation leaves open the question whether the hnite 
CO- dimensional manifold of radial hnite energy data scattering to unstable steady states is 
closed in the energy topology. The answer to this question seems to be nontrivial and will 
require further understanding of the global dynamics. For example, consider an unstable 
excited state (0,0) G S. It is not hard to show that there is a radial solution lt{t) which 
converges to (0, 0) exponentially as f —)■ —oo, i.e., lt{t) is on the unstable manifold of (0, 0), 
and hence £{lt{t)) = By the channel of energy property established below, lt{t) 

will emit a nontrivial amount of energy to large distances as t —?■ -|-oo and subsequently 
scatter to a steady state of strictly less energy, say (0, 0). However, there is a possibility that 
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(0,0) is also an excited state. In that case, denote by the manifold of data scattering to 

(0, 0) as t ^ +CX 0 , we see (t) G for all t, but it (t) (0, 0) in x as t ^ —oo and 

clearly (0, 0) ^ Consequently, in such a situation would not be closed. Admittedly, 

such behavior should be non-generic (due to the fact that 0 might be expected to be the 
ground state) and perhaps impossible for a generic choice of V. We plan to address this 
question in future work. 

Another interesting question is if this description of global dynamics can be achieved 
without the radial assumption. This question seems to be very challenging. Recall that 
in the radial case, due to the channel of energy property, we only need to consider the 
dynamics outside some well-chosen light cone, where the dynamics is relatively simple. In 
contrast, in the non-radial case, where only less effective channel of energy inequality is 
available, one must deal directly with the complicated dynamics in a finite region. In this 
case the only other global tool is the virial type identities. However the presence of the 
spatial inhomogeneity V seems to render such identities ineffective. In particular, we do not 
know if the only compact solutions are steady states. Recall that lt{t) is called compact 
if {^(t) : t G M} is pre-compact in x L^. This is in sharp contrast with the energy 
critical focusing wave equation case, where one knows exactly what these compact solutions 
are (modulo some non-degeneracy condition on steady states). Hence, a full characterization 
of compact solutions seems to be a natural first step. 

Our paper is organized as follows. In Section [2l we study steady states to equation fll.ll) 
and show in particular the existence of stable excited steady states; in Section [3] we construct 
the local center-stable manifold. The novelty of his construction lies with the fact that it 
is carried out near any solution which scatters to a given unstable steady state, without, 
however, being necessarily close to the steady state in the energy topology, in Section 0] 
we recall some results on the well-known profile decompositions and channel of energy in¬ 
equalities, adapted to equation fll.ll) : in Section [5] we prove our main result Theorem 11.31 
Appendix A contains some elliptic estimates for the steady states; Appendix B proves an 
endpoint Strichartz estimate for the inhomogeneous wave equation in the radial case. 


2. Steady state solutions 


In this section we prove some results about the steady states that are relevant for the global 
dynamics. We first give necessary and sufficient conditions for the existence of nontrivial 
ground state. Recall that such a state is the global minimizer of the energy functional 

J(0) := f 


iv(?r 


2^6 


dx. 


( 2 . 1 ) 


Lemma 2.1. Consider J as a functional defined in If the operator —A — V has 

negative eigenvalues then there exists a global minimizer Q > 0 with J{Q) < 0. If —A — V 
has no negative eigenvalues, then the only steady state solution u G to eguation / li.i|) 

is u = t). 


Remark 2.1. The proof of this lemma is a simple application of variational arguments and 
the strong maximum principle, we omit the standard details. 
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In the case that —A — V has no negative eigenvalues and assuming that we only consider 
radial solutions, then from the results in [16] we know that all radial hnite energy solutions 
to equation fll.ip scatter to the trivial steady state. In what follows we therefore assume 
that —A — V has some negative eigenvalues, so that we have nontrivial global minimizers Q 
and —Q. We call Q and —Q ground states, and call other steady solutions excited states. 


The next result shows the uniqueness of ground states. Note that we do not need radial 
symmetry here. 

Lemma 2.2. There is at most one nontrivial nonnegative steady state in to eguation 

anp. 

Proof. Suppose Qi and Q 2 are two nontrivial nonnegative steady solutions to equation fll.ip 
in Then 

-AQi - l/Qi + Q? = 0, 

~ A(52 ~ V Q2 + Q2 — 0. 

By standard elliptic estimates, we have Qi, Q 2 G P < 00 . Moreover we 

have the following decay estimates 

|Qi(x)| + \Q2{x)\ < ^ , X e (2.2) 


1 + X 


The above claim on the regularity and decay holds for any steady state in and follows 
from more or less standard elliptic techniques. For the sake of completeness, we provide a 
proof in the Appendix A. By Strong Maximum Principle, we see that Qi, Q2 > 0. Denote 
the open set 

D := {x : Qi(x) > Q 2 {x)}, 

we have 


ft (-A<3i - VQi + <3?) - Qi (-AQ 2 - VQ 2 + Q^) dx = 0. 


(2.3) 


By the regularity and decay properties of Qi, Q 2 we can integrate by parts in equation fl2.3p . 
noting that Qi = Q 2 on dVL, we see that 


Note that 


and 


f Qi-^{Q 2 - Qi) da + f QiQ 2 {Qt - Qt) dx = 0. 
'an on Jq 


{Q2 - Qi) > 0 on cAl, 
on 


(2.4) 


Q1Q2 (Qt - Q2) >0 ill 

Thus equation fl2.4p can hold only if D = 0. Thus Qi < Q 2 . Similarly Q 2 < Qi. Therefore 

Qi = Q2. □ 

Naively one might expect excited states to be unstable, since they change sign. However 
in general this may not be the case, as seen from the following theorem. 












Theorem 2.1. There exists an open set O G V such that for any V G O, there exists an 
exicted state 0 to equation which is stable. 


Proof. The proof is based on simple perturbation arguments, once a good large potential is 
chosen. We can construct an excited state near a good “prohle”, for which the linearized 
operator is explicit and stable, with a well-chosen potential. Then we can conclude that the 
linearization near the excited state is also stable. 

Step 1: choice ofV. Denote 



as the unique (up to scaling and sign change) radial solution to 

- Au = u^. (2.5) 

Let us take 

l/l(x) = 2W^ (2.6) 

and positive A sufficiently large to be chosen below. Set 

Du(a;) = AVi(Ax). (2.7) 

It is easy to check that W solves 

—Am — Viu -I- M® = 0, 

and that W\{x) = X^W{\x) solves 

— Am — Vi\u + = 0. 


We choose V := Vi + Wa- 


Step 2: construction of a stable excited state. Consider the following elliptic equation 


- A0 - 1/0 + = 0. (2.8) 

Our goal is to show if A is sufficiently large then we can construct a steady state 0 of the 
form 

^ = W-Wx + T], (2.9) 

with some small r]. The equation for p is 

-Ap + {3W^ + 3W^ - 20WW^ - 20WxW^ + 30W^W^)p + N{p, A) = /a, (2.10) 

where the nonlinear term is 

N{p, A) = 10(W - ITa) V + 10(W - PLa) V + 5(1D - Wx)p^ + (2.11) 

and the nonhomogeneous term is 

fx = -(PL - ILa)"^ + PL'^ - PP^a' - PAPLa + PLPAa- (2.12) 


If A is sufficiently large, fx will be small in appropriate function spaces and we can use a 
perturbation argument to solve for p. A key ingredient is the following standard uniform 
estimate in A on the linear part: 
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Claim 2.1. For sufficiently large A, the operator 

La := -A + 3PF^ + 3W^ - 20WW^ - 20PFa1^^ + SOPF^PFa^ : H\R^) 


is invertible and we have the following estimate on the norm of the inverse operator L^ 


-1 


'-ii 


IIh-1(R3)-i>H1(R3) 


<C, 


(2,13) 


where C is some absolute constant. 


Proof of Claim f^TZl It is easy to verify that —A : Lf^(R^) —)■ is invertible, and that 

La is a compact perturbation of —A. Thus to prove that La is invertible we only need to 
show its kernel is trivial. This follows directly from the following bound for large A 

e (2.14) 

where the inner product is with respect to the and H~^ pairing. The proof of fl2.14l) is 
an easy consequence of integration by parts argument and Holder’s inequality, once we note 
that 

lim IllTlT^ + ^ W^W^\\L3/2r^3) = 0. 

Suppose L\u = /. Then from the bound fl2.14p and Holder’s inequality we infer that 

II Vm||l2(ir3 ) < 2||/||^-i, for sufficiently large A. (2-15) 

In view of the preceding, the bound 02.131) follows immediately. 

Using the uniform bound on we can easily solve for r]. 

Claim 2.2. For any e > 0, if A is sufficiently large, then there exists a solution rj G Lr^nL®(M^) 
to equation O2.10p . in the sense of distributions, with 

||hlU6(K3) < Ce, (2.16) 

where C is an absolute constant. 


Proof. Take small e > 0. We can take A sufficiently large, so that ||/a||l6/5 < e and 02.131) 
holds. We reformulate equation O2.10p in the following way 

V = Lf^fx-Lf^Nir^,X). (2.17) 

Note that embeds continuously into thus the right hand side makes sense. 

Now one can check that 

ioVA-iyivo.A) 

is a contraction mapping in i? 2 e C L®(M^), if we choose e small enough. Thus equation 02.171) 
and consequently equation 02.101) has a unique solution r], with ||?7||L6(ffi3) < 2e. Clearly, this 
7] satishes the requirement of Claim 2.2. 

By looking at the L® norm of the positive and negative parts of cf in 02.Op . it follows that 
the steady state W — Wx + rj changes sign if we choose e sufficiently small, and thus is an 
excited state. Moreover the linearized operator around this excited state is 

-A - 21F^ - 21Fa^ + 5{W - Wa + r/)^ 
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If we choose e small enough one can show this operator is nonnegative and has no negative 
eigenvalues nor zero eigenvalues/resonance. A standard local perturbation analysis implies 
that the excited state W — Wx + rj is stable. 

Step 3: stability with respect to V = Vi + Vix. We now show that our construction is stable 
with respect to small perturbations of the potential V in Y. This is more or less clear from 
the existence proof. Below we just outline some key points. We now formulate 

Claim 2.3. Let V := Vi + Wa be dehned as above. Assume A sufficiently large and 6 
sufficiently small, then for any radial potential V &Y satisfying 

||C-C||y <5, (2.18) 

there exists a stable excited state to 

- A0 - 1/0 + = 0 in (2.19) 

Proof of Claim \2.3i For any V eY with 

IlC-Clly <5, 

for some 6 sufficiently small. We have 

||(_A-C)-(-A-C)||^,^^_, <C<5. 

Hence if we choose A sufficiently large, 6 sufficiently small, we can repeat Step 2 with V 
replaced by V. The resulting Lx will still satisfy the uniform bound in Claim 12.11 fx with 
extra terms iV — V)W and {V — V)Wx, and N{ri, A) can still be controlled in exactly the 
same way, as long as 6 is chosen sufficiently small. One can then use the Contraction 
Mapping Theorem to hnish the proof. We omit the routine details. This hnishes the proof 
of Theorem 12.11 

□ 


The above theorem proves the existence of stable excited states, on the other hand it is 
clear that there are unstable excited states. For example, suppose that the operator —A — V 
has negative eigenvalues (so that there are nontrivial ground states), then the excited state 
0 = 0 is unstable. A perhaps more interesting fact is that “newly bifurcated” excited states 
are unstable. 


Lemma 2.3. Let a G [0, oo), V E C“(R^) be nonnegative and not identically zero. Suppose 
that for a > ai the principal eigenvalue Ai(q!) of the operator —A — aV is negative. Assume 
further that 

(^ 2 , 0) E («!, oo) X Fr^(R^) 
is a bifurcation point for the eguation 

- A0-aI/0 + 0® = 0 , ( 2 . 20 ) 


in the sense that for any e > 0 there is a nontrivial radial steady state (a, 0) to 
^KM) with 


a - a 2 \ + ||0|lifi(R3) < e. 


Then for e sufficiently small, the steady state {a, 0) satisfying H2.21]) is unstable. 


eguation 

( 2 . 21 ) 
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Proof. Choosing e sufficiently small we have a > ai. Since V is nonnegative, we have that 

'^i(ci) ^ -^i(cii) ^ 0, 

thus there exist 5 > 0 and rjj e such that 

f {\V^l^\^-aV^l^^)dx< j ( 2 . 22 ) 

J J 

Note the linearized operator around (f is 

- A - aC + 5(p\ (2.23) 


If e is sufficiently small, by Holder’s inequality, we obtain 

f {\V^|J\^-aV^|J^ + 5(|)'^^f^)dx<-i6-Ce^ 


hi 


< 0 . 


Thus the linearized operator has at least one negative eigenvalue and consequently the small 
excited states are unstable. □ 


Remark 2.2. For a description of bifurcations, see Appendix A of [16]. The above arguments 
also imply that if —A — V has negative eigenvalues, then all sufficiently small excited states 
in are unstable. 

In Theorem 6.1 of [IB], we showed that there exists a dense open set Hi C y such that for 
any 1/ G Hi, there exists only hnitely many radial steady states to equation fll.ip . all of which 
are hyperbolic when the linearized operator is restricted to the space of radial functions. As 
is mentioned in the introduction, this leaves open the possibility that some radial steady 
states may still have zero eigenvalues or a zero resonance without the radial assumption. As 
the property of having zero eigenvalues or a zero resonance is non-generic, we can expect to 
eliminate such behavior by removing from Hi a closed set while making the remaining set 
still open and dense. More precisely we have the following result. 

Lemma 2.4. Let Hi be a dense open subset of Y such that for any 1/ G Hi there are only 
finitely many radial steady states to equation all of which are hyperbolic if the linearized 

operator is restricted to radial functions. Then there exists a dense open set H C Hi, such 
that for any H G H, all radial steady states are hyperbolic without radial symmetry. 


Remark. We hx this choice of H below. 

Proof. Take any V G Hi, and suppose that 0i, (j) 2 , ■ ■ ■ ,4>n are the radial steady states corre¬ 
sponding to V. Since the linearized operator around any (fi is hyperbolic in the space of radial 
functions, standard perturbation arguments imply that any sufficiently small perturbation 
of y in H will not change the number of radial steady states, and each radial steady state (pi 
depends smoothly on the perturbation. We will show that we can hnd a specihc perturbation 
V oi V with arbitrarily small norm in H, (in particular V G Hi), such that the perturbed 
steady state piiV) becomes hyperbolic even without restricting to radial functions. Since 
being hyperbolic is an open property, we can then make repeated small perturbations to the 
potential until all the steady states become hyperbolic even in the nonradial function space. 
Then it is clear the subset H C Hi with the property that for any 1/ G H, all the radial 
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steady states are hyperbolic without restriction to radial functions, is dense and open in Y. 
Below we describe the perturbation in detail. For e > 0 sufficiently small, set 

V-=V + e(p{. (2.24) 

Suppose that the perturbed steady state becomes 

(j)i := (j)i + ei/j. (2.25) 


Then the equation for ip is 

— Apj — Vp: + = PI + e piip — N{'ll:, e), in (2.26) 


where 

N{Tjj, e) = lOe + e^ip^. (2.27) 


Noting that 

-A01 - Vpi + 501 • 01 = 40f, 

and the fact that the linearized operator —A — V + 50f is invertible from ^ H~^r\L‘^ 

when restricted to radial functions (which follows from hyperbolicity in radial functions), we 
can rewrite equation fl2.26p as 


0 = ^ + (-A - 4^ + 50t)-i(e 0t 0 - A(0, e)). (2.28) 

If we take e sufficiently small we can assume 17 G hli, and moreover we can use standard 
perturbation arguments to show that 


Then the linearized operator around the perturbed steady state 0i + e 0 becomes 


- A-1/ 

The key point for us is that 




0 + 0 3 (e^ 

1 L7nL°°^ 


5(1 + 1) -e 


01 > 501 + 4e 01, 


(2.29) 


(2.30) 


(2.31) 


hence we have gained a positive factor which will eliminate the zero eigenvalues/zero reso¬ 
nance. The proof is finished with the following claim and the Min-Max principle for eigen¬ 
values. □ 


Claim 2.4. Let V and 0i be given as above. Suppose that the linearized operator 

-A - 1/ + 50t 

has k > 0 negative eigenvalues with corresponding eigenfunctions pi ,..., Pk, and possibly 
also zero eigenvalues or zero resonance. Then for e > 0 sufficiently small, and any / G 
with 

/ fp^dx = ---= fpkdx = 0, 

7r3 7r3 








we have 
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/ (-y + 5<Pt + + \Vf\^dx > ce \\f\\%„ 

JK3 

for some fixed c > 0. 

Proof. Consider the functional 

A(/):= [ {-V + 5(ft)f+\^ffdx 

JR3 

on the space 

X:=lfeH^: I fpidx = ---=l fpkdx = 0}. 

Suppose 


(2.32) 


(2.33) 


(2.34) 


feH^-. (-A-C + 50t)/ = o} =span{Zi,...,Z4 


Using the same arguments as in Proposition 3.6 in [12], we can find linearly independent 
Ei,...,Eme C^, with 


Vz = 1... fc, Vj = 1... m, y piEj = 0, Vz, j = 1... m, y (f^Ei Zj = 6ij, 
such that for any / with 

[ / Pi = 0, [ f Ej(j){ = 0, Vz = 1... /c, j = 


one has 


For any / G X, we can decompose 


A(/) > c 


ifi- 


(2.35) 

(2.36) 

(2.37) 


f = '^CiEi + g, with j gYi = ■■■ = j gYk = 0, j gEi(f{ = ■■■ = j gE^(j)l = 0 . 


We distinguish two cases. 


(2.38) 


Case 1: ^ \ci\ > for some small > 0. Then the claimed bound fl2.32p follows 

i=l 

from 


/ „ m m 

/ (^ CiEif (f){dx>e^ c- > 
d i=i i=i 


2 

pl’ 


(2.39) 


for some small 82 > 0, where in the third inequality we have used the equivalence of norms 
in a hnite dimensional space and the linear independence of Ei. 

m 

Case 2: \c\ ^ ll/llpi- Then the claimed bound fl2.32p follows directly from the bound 

dMZI). □ 
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3. Construction of the local center-stable manifold 


In [TB] we showed that there exists a dense open set Q G Y, such that for any potential 
1/ G there are only hnitely many radial steady states and all radial hnite energy solutions 
scatter to one of the steady states. Furthermore, the linearized operator —A — V + 50^ has 
no zero eigenvalues nor zero resonance for any radial steady state (j). Our goal in this section 
is to study local dynamics around a solution which scatters to an unstable excited state. 
Recall 

S{t) = {So{t),S,{t)), teR 

denotes the solution flow. Thus given (mo,Mi) G x {t)(uo,Ui) is the solution 

to equation fll.ip with initial data (mo,Mi). To prove our main result, let us hrst state the 
Strichartz estimate for the linear wave equation in dimension 3. 


Theorem 3.1. Let I be a time interval and let v : I xR? ^ R he a finite energy solution to 
the wave equation 

{du -A)v = F 

with initial data {v(to),dtv(to)) = {f,g) for some to G I. Then we have the estimates 

C(9, r) (|I(/,S)|| hixL2 + LiLi(/xR3)) , (3.1) 

where 2 < q < oo and 2 < r < oo satisfy the scaling condition 


1 + 3 _ 1 
q r 2 

and the wave admissible condition 

1 1 1 

q r 2 

Moreover, if f,g and F are radial functions, li3.1\) holds true when {q,r) = (2, oo). 


(3.2) 

(3.3) 


Later, we will call a pair {q,r),2 < q,r < oo admissible if it satishes fl3.2p . 03.31) . 

See [13] for the proof of this theorem in the nonradial and non-endpoint case. The for¬ 
bidden endpoint {q,r) = {2,oo) was found in |20|, where it was also proved that for the 
homogeneous equation [F = 0), (2, cx)) becomes “admissible” if the initial data is of the 
form {v{to), dtv{to)) = (0, g) with g radially symmetric. For the sake of completeness, in Ap¬ 
pendix B, we provide the proof of the endpoint Strichartz estimate for the inhomogeneous 
equation with general radial data. 


For applications below, we need the following Strichartz estimates for solutions to linear 
wave equation with potential. 

Lemma 3.1. Take V E Y such that the operator —A — V has no zero eigenvalues or zero 
resonance. Denote as the projection operator to the continuous spectrum of —A — V. 

Denote _ 

a; := ^yP^{-A - V). (3.4) 

Let I be a time interval with to E F Then for any (/, g) E H^xLf{R^) and F E LlL‘^{IxR^), 
the solution {t) to the equation 

dtt'y + = P-^F, {t, x) E I xR^, 


(3.5) 
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with ^(to) = satisfies 

\\i.int)\\c°{H^xL^) + I|7||l?L5(/xR3) < (||(/)5')IIh1xL2 + II-^IIlIlk/xrs)) , (3.6) 

where 2 < q < oo and 2 < r < oo satisfy the conditions and Moreover, if f,g 

and F are radial functions, H3.6\) holds true when {q,r) = (2, cxo). 


Remark. Our proof of Strichartz estimates fl3.6p relies crucially on U’ boundedness of wave 
operators with minimal decay conditions on V, obtained by Beceanu [1]. For earlier impor¬ 
tant work on boundedness of wave operators, see Yajima [30] and references therein. 

Proof. The energy estimate on ^ follows from standard integration by parts argument and 
inequality fld.lUp below. We will therefore concentrate only on the Strichartz estimates. 
Recall the definition of the forward wave operator 

W+ := s - lim e^h-A-v)^-iti-A)_ 

t—>00 

It is well known that for V G Y, W+ is linearly isomorphic from to see for 

example [2]. In |1] Beceanu obtained among other things an important structural theorem 
for hF+, which implies that hF+ is also bounded in U’ spaces with 1 < p < 00 : 

WMlv<\Mlv VpeL^nLF (3.8) 


Hence 1F+ can be naturally extended as a bounded operator in for 1 < p < 00 and in 
which is the completion of Lf nL°° in L°°. An important fact of W+ is the intertwining 
property 

= WMmWl (3.9) 

which holds for p G and also for more general <p by limiting arguments provided that 
one can obtain suitable bounds. Here Wf is the adjoint operator of W+. Note that {t) 
admits the following representation 


7(t) = cos (cut) P^f -I- 


sincut 


P"g+ / 
^ Jo 


sincu(t — s) | 


cu 


= W+cos{\V\t)Wff+ + W. 

= / + // + ///. 


P^F{s) ds 

* sin |V|(t — s) 


'0 


WfP^Pis) ds 


The Strichartz estimates for part H and HI then follow directly from the boundedness 
fl3.8p of IY+ and the corresponding Strichartz estimates for free radiations. For the endpoint 
(g, r) = (2, 00 ), we only need to note in addition that leaves the space of radial functions 
invariant, since V is radial. It remains to consider I. Firstly we claim 

Claim 3.1. The operator cu which is initially defined in H^, satishes 

l|V»^|li. < j e 

and can be extended naturally as a linear isomorphism from to P^L^. 


(3.10) 
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Let us assume this claim momentarily. Then the Strichartz estimate for part I is easy to 
prove. By Claim IXTI we can write P-^f = uj~^f for some / G P-^L^. Thus 

W+cos {\V\t)Wlf = W+cos {\V\t)Wluj-^f 
= lT+cos(|V|f)|V|-ifT:/, 

where we have used the fact that on P^L^, which follows from the 

intertwining property fl3.9p by suitable limiting arguments with the help of bounds fl3.10p and 
fl3.8p . Note that W^f G consequently |V|“^1T+/ G H^. Hence the Strichartz estimates 
for part I follow straightforwardly from the corresponding estimates for free radiations and 
the bound fl3.8p . □ 

Now we give a brief proof of Claim ITTI The second part of the inequality fl3.10p is an easy 
consequence of the fact that u is self adjoint and = P^{—A — H), and an integration by 
parts argument. The first part of the inequality follows from the assumption that —A — V 
has no zero eigenvalues or zero resonance. Hence oo can be extended as a bounded operator 
from P^H^ to P^L^. Moreover u has closed range due to the bound fl3.10l) . Since —A — V 
has no zero eigenvalues, we conclude that the range of A — H) is dense in P^L^. Thus 
the range of uj which is bigger than the range of A — H) is also dense in P^P^. Com¬ 
bining these two facts, we see that uj is indeed a linear isomorphism from P^H^ to P^L^ 
and the claim is proved. 


Our main goal in this section is to prove the following result. 

Theorem 3.2. Let Q be a dense open subset ofY such that equation M.l\) has only finitely 
many radial steady states, all of which are hyperbolic. Suppose V ^ Vt <zY. Suppose Tf (t) 
is a radial finite energy solution to equation M.l\) which scatters to an unstable steady state 
(0,0). Let 

-kl<-kl<---<-kl<0 (3.11) 

be the negative eigenvalues of —A — V + 50^^ restricted to radial functions (counted with 
multiplicity) with normalized eigenfunctions pi, P 2 , ■ ■ ■ ,Pn, respectively. Decompose 

-^rad ^ -^rad(®^) = © ^u, (3.12) 

where 


Xs = \ {uo,ui) G P^ad X -^Ld(®^) : {kjUQ + mi,Pj)l 2 = 0, for all 1 < j < n| , (3.13) 


and 

Xu = span {{pj, kjPj), 1 <j <n} . (3.14) 

Then there exist Cq > 0, T sufficiently large, a ball Pep((0,0)) C x L^^j(M^), and a 


HLd X Lrndy 


(3.15) 


smooth mapping 

T: f^(r) + (P,a((0,0))nA,) ^ 

satisfying T(L^(T)) = iJiT), with the following property. Let M be the graph o/T and set 
M. = ^(—r)A4. Then any solution to equation U.l\) with initial data {uo,ui) G Xi scatters 
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to (</>, 0). Moreover, there is an ei with 0 < ei < eg, such that if a solution it it) with initial 
data (mo,Mi) G -Be^(t/(0)) C x satisfies 

\\lt{t) — iJ (i)||i 7 ixL 2 < ei for all t >0, (3.16) 

then {uq, Ui) G Ai. 

Proof. By assumption, there exists free radial radiation such that 

hm ||ft(() - (^,0) - h'-(i)llA.xy = O- (3-17) 

t^OO 

We divide our construction of the center-stable manifold into a series of steps: 

Step 0: L® decay for free waves. We observe that for any finite energy free radiation it^, we 
have 

||'*^^(^)I1l 6 0 as t —)■ cxo. (3.18) 

For smooth t^^(O) with supp^^(O) <1 Br, we have 

C 

\u^{t,x)\ < -jXt-R<\x\<t+R, for t>R. 

Then fl3.18p follows by direct calculation. For general initial data, fl3.18p follows from approx¬ 
imations by compactly supported smooth functions and the uniform bound ||M^(t)||i6(R3) < 

C||^'^(0)||p^lxi2(K3). 

Step 1: space-time estimates for U — fi. Denote h{t, x) = U{t, x) — 4>{x), then h satishes the 
equation 

htt — Ah — V{x)h + Icjfih -|- A(0, h) = 0, (3.19) 

where 

Ar(0, h) = {(!) + hf -(ft - bcjth. 

Since iJ G T“([0, cxo), x L^) and U G x M^) for any finite interval /, by equation 

fll.ip and Strichartz estimate fl3.ll) we see U G x M^). By standard elliptic estimates, 

we know that 0 G C^(M^)(see also Appendix A), hence h G x R^) for any finite time 

interval I. In what follows, we will show that 

ll^llL^Lg°([0,oo)xM3) < OO- 

Recall that pi,..., are the n radial Lfi normalized orthogonal eigenfunctions of the oper¬ 
ator 

C^ = -A-V + 50 ^ 

corresponding to the eigenvalues (counting multiplicity) —k\ < ... — < 0, respectively. 

From Agmon’s estimate [1], we know these eigenfunctions decay exponentially. Writing 

h = Ai(t)pi \n{t)Pn + 7) 

with 7 T Pi for i = 1, • • • , n, and plugging this into equation fl3.19l) we obtain 

n 

^(Ai(f) - kiXi{t))pi -h 7 -h £^7 = N{(j), h). 

i=l 


(3.20) 
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Denote by P* the projection operator onto the th eigenfnnction and by the projection 
operator onto the continnons spectrnm restricted to radial fnnctions0, i.e., 


Pi = Pi® Pi, 


P^ = I 


i=l 


Pi ® Pi 


Applying the projection operators Pj and P-*- to eqnation (I3.20p . we derive the following 
eqnations for Xi{t) and 'j(t,x): 


Xi{t) - k^Xiit) = PiN{(j), h) := i = l,...n 


7 + 01^7 = P-^N{(j), h) := Ac, 


a; := ./P^£, 


(3.21) 


Note that the steady state 0 decays at the rate as |x| —)■ cxo, hence the potential 


in 


the operator which is —V + 50^, decays like Q( (-]^_)_|,^|^Ln{/ 3 , 4 } )• This decay rate is better 
than the critical rate 0{j^) as |x| —)• oo (in fact —V + 50^^ G Y). Hence we can apply 
the resnlt of Lemma 13.11 and conclude that Strichartz estimates as in Lemma 13.11 hold for 
solutions of the equation 

Vtt + = F, (3.22) 

with F radial and satisfying the compatibility condition P-*“P = P. 

From fl3.17p and fl3.18p . we know that 

lim \\h{t, x)\\l?°lU[t,oc)xr^) = 0 

t^OO I J.VI 

Also using the fact that pi decay exponentially, we have 


\Ht)\ = \{Pi\h)\ < ||pi||^6||h(t,a;)||i6(K3) ^0 


as t ^ oo. 


Let r(t) be the solution operator to the equation Vu + oo'^v = 0, i.e.. 


T{t - to)(7(^o),7(^o)) = cos(a;(t - to))7(^o) + -sin(a;(t - to))7(^o) 

00 


We claim: 


Claim 3.2. Given any e 1, we have 

||r(t - t)(7(t),7(t))|1^2^^([7.oc))xe 3) < e. 

for sufficiently large T > 0. 

We postpone the proof of Claim 13.21 to the end of the proof of Theorem 13.21 
Hence given a small positive number e 1, which will be chosen later, we can pick a large 
time T = T{e,U), such that 

\\h\\ l^lU[t,oo)xr 3) Y e (3.23) 

ll^i(^)llLj°“([T,oo)) < e (3.24) 

||r(t - T)(7(T),7(T))||^2£,oo([-roo)xK3) < e. (3.25) 

'^Note that can be written via Stone’s formula as integral of resolvant of hence it is invariant for 
radial functions. In particular, does not involve any non-radial eigenfunctions. 
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From the equation for Aj(t) in fl3.2ip . we conclude that for t >T 


Xi{t) =cosh.{ki{t - T))Xi{T) + — sinh(/ci(f - T))Ai(T) 

r^i 


ki Jj 


sinh(fcj(t — s))Np^{s) ds 

i\y{ JT 


ek,iT-s)N^^^s) ds 




where 'R-{t) denotes a term that remains bounded for bounded Np.{s). By (13 .24^ . we obtain 
the following stability condition 


A,(r) = -kMT) - / e^^^'^-^^Np^{s)ds. 


(3.26) 


Under this condition we can rewrite equation fl3.21l) as the following integral equation 

1 /■°° .1 1 


\.{t) = 


A,(r) + — / 


i JT 


2 k 


e-’^^\^-^\NpXs)ds, 


i JT 


7 (t) = cos(a;(t - T))'y{T) H— sin(a;(t - T))'j{T) + 


sin(a;(t — s)) 


(3.27) 


cn 


u 


Nc{s)ds. 


For any time T > T, we dehne 

n 

||(Al, . . . , An, 7)||x([r,f)) • ~ ll'^*(^)llL2([7’^'f)) + WxW L^L'^{[T,f)xR3)- 

i=l 

With the help of Strichartz estimates from Lemma [3.11 by estimating fl3.27p . we get that 
ll-^i(^)llL2([r,T)) —^1 + ll-^PillLi([r,T)) IIl7=([t,oo))) ’ (3.28) 

\\x\\L'^L^([T,f)xR3) —^2 (^||r(t “ ^)(7(^), 7(^))lli:,2ioo('[2^f)xR3) + \\^c\\i,iLl{lT,f)xR^)) ' (3-29) 

Here that constant Ci depends on the and integrals of and the constant C 2 
depends only on the constants in the Strichartz estimates. 

In fl3.29p . instead of estimating initial data ( 7 (T), 7 (T)) in x L^(R^) which may not 
be small, we estimate its free evolution in LfL'^{[T,T) x M^). Consequently, we obtain 
smallness because of fl3.25p . 

Using the fact that 

Np^ = (p,|iV(0, h)), Np = Y^ Np^p,, N, = N-Np 


and the exponential decay of pi, we have 

\\N,M Li([T,r))> ll^c|lLiL2([T,r)xR3) < llLiL2([T^r)xR3)- 

Recall that 

N{(p, h) = + 50h^ + h^. 


(3.30) 
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Hence, by Holder inequalities, we have 

\\LlLl{[T,f)xR3) < 11‘^lli® ll^llL2igo([T,f)xR3)’ 

11*^ ^ \\LlLl{[T,f)xR3) < ll^llL2Lgo([T,f)xR3)ll^lln?°'^E(['r-T^)x®®)’ 

lll,jL2([T,f)xR3) — W^W^xW^W L'^L'^([T,f)xR3)\\^\\L^L^(^[T,f)xm3)^ 

11^ llLiL2([T,f)xR3) ^ \\^\\L^L'^([T,f)xR3)\\^\\L^LU[T,f)xR3)^ 

Consequently 

11^(05 ^)llLiL2([T,f)xR3) - ^\\H'L‘2L^{[T,f)xM.S)^ (3.31) 

and 




< CE ll•#-ll&il'>lli^.i|(l^.=o,xR., £ C:e\ (3,32) 

i=2 

Using (13.2311 and 

ll^llL2Loo([r,T)xR3) r-^ II(Ai,--- )A„, 7 )II 

with constant depending on ||pdlL°°, we can combine estimates (13.281) . (13.291) with (13.241) . 
(Km and (I330D - (I332D to get 


II (Ai, • • • , An, 7 ) llx([r,T)) — + II (Ai, • ■ ■ ) An, 7 ) + e I, 

Here K is some constant depending only on the constants in Strichartz inequalities for 
equation (I3.22p and ||0 ||l 6 and ||pj||Lgo. Since this estimate is true for all T > T, we can 
choose e <C 1, which can be achieved by taking T sufficiently large, such that 

{AK^ + l)e < 1. 


By a continuity argument we then obtain that 

||(Ai, • • • , An,7)||v([r,cx))) < 2 Ke. 

which implies that ||/i||L 22 ,go([T,oo)xR 3 ) Using interpolation between the and 

norms, we can also obtain 

II^IIl?L 5 ([t,oo)xr 3 ) < e, for any admissible pair (g,r),g > 2. (3.33) 

In particular, we infer that h G LjL^°([0, 00 ) x M^). 

Step 2: construction of the center-stable manifold near a solution U. Given a radial finite 
energy solution U to (II.Ih satisfying (I3.17p . we consider another radial finite energy solution 
u, with II (T) — it (T)||^ 1 xl 2 (r 3 ) small for a fixed large time T from Step 1 (we may need to 

take T large to close the estimates below, which of course can be done). We write u = U 
then 7 satishes the equation 

r]u - A 7 - V{x)r] + {U + pf - = 0, (t, x) G (T, 00 ). 
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Plugging in 17 = 0 + /i, we can further write the equation as 

% + £077 + iV(0, h,r 7 ) = 0, {t,x) e {T,oo), (3.34) 

with 

N{(j), h,r]) = {(j) + h + vif - (0 + hf - 5(j)‘^r] 

We note that N still contains terms linear in rj. However, a closer inspection shows that the 
coefficients of the linear terms in rj decay in both space and time, and can be made small if 
we choose T sufficiently large. First write 

T] = Xi{t)pi H-h \n{t)pn + 7, 7 -L A 


for 7 = 1, • • • , 77 . We use similar arguments as in step 1 to obtain a solution 77 which stays 
small for all positive times with given (Ai(T),-- - ,A„(T)) and We can obtain 

equations for A ,,7 similar to fl3.2ip . Since we seek a forward solution which grows at most 
polynomially, we obtain a similar necessary and sufficient stability condition as fl3.26p 


i(T) = -kMT) - j 

Using equations 03.341) and 03.351) we arrive at the system of equations for Aj and 7 , 


\i{t) = e 




HT) + 


2A 




1 


H JT 

1 • 1 /■* 

7 (f) = cos(a;(t — T)) 7 (T) H—sin(ci;(f — T)) 7 (T) H— / sm(u{t — s))Nc{s) ds. 

U (jJ Jt 


Dehne 


||(Ai,..., An,7)||x ^ \\Xi{t)\\L^riLl{[T,oo)) + ll7|lL7>°h'inL2L^([r,oo)x 


i=l 


Estimating system 03.36p . we obtain that 

||Ai(f)||L°°nL2([r,oo)) ^|Ai(£)| + ||^Pi||Lj([T,oo)) ^ l^i(£)l + II^IIl1l2([T,oo)x 
ll7llL7>°h'inL2L7=([T,oo)xR3) ^ll(7(£),7(£))lli7-ixL2 + II^IIl1l2([T,oo)xR3)- 
Recalling |iV| < \(l)^-^h^r]\ + 'Ek> 2 ,i+j+k =5 we have 

LjLl{[T,oo)xm.3) < ll^llie ||h||j;^22,oo(['r^oo)xR3)||77||j;^2^oo([7-_oo)xR3), 

\\<P‘^h‘^V\\LlLl{lT,oo)xK.3) < ||0||L6|l^llA°°L6([T,oo)xR3)||h||2,22,oc([-r^^)x 

\\4>h^v\\LlLl{[T,oo)xm.3) < \\4>\\L^\\h\\L^LU[T,oo)xR3)\\HL'^L'^{[T,oo)x 


L2Lgc([r,oo)xI 

L2l-([T,oo)xI 


II^^7||l1L2([T,oo)xR 3) < ll^lli7>°L6([T,oo)xR3)II^IU2Lg°([T,oo)xR3)||7llL2Lg°([T,cx))xR3)- 


Using 03.331) . we get that 

4 




i=i 


< 


L2L-([r,oo)xIt 


LiL2([T,oo)xR3) 


(3.35) 




(3.36) 


(3.37) 

(3.38) 

(3.39) 


(3.40) 
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The higher order terms are easier to estimate. We can always place h in whence 


Since 


k>2^i-\-j-\-k=b 


L'fmr\^L^{[T,oo)xl 


5 


Lf^ H^nLlL^({T,oo)xl 


(3.41) 


LiLi([r,oo)xR3) ^-2 

< ||(Ai,--- , A„, 7 )||x([t,oo)), we can combine the preceding 


estimates and get that 

ll(Ai, ■ ■ ■, a„,7)IU([t,„)) < Af ( E + ll(A(r).4(r))llt,.xi> 


. 2 = 1 


+ -^'^<a||(Ai, • • • , An, 7)I|2 s:([t,oo)) + K Ell7i.'"4«.7)llW,co),. 


k=2 


where ii' > 1 is a constant only depending on the constants in the Strichartz estimates for 
equation fl3.22p and || 0 ||l 6 (r 3 ) and ||pi||i,j°. This inequality implies that if we take e = eo 
sufficiently small (which can be achieved by choosing T suitably large), and 6 < Cq with 


5;|A,(r)| + ||(7(T),yr))||ft,,i, <i, 


(3.42) 


2=1 


such that Keo < 1 and K'^S < then the map dehned by the right-hand side of system 
fl3.36p takes a ball B2 ks{0 ) into itself. Moreover, we can check by the same argument that 
this map is in fact a contraction. Thus for any given small (Ai(T), • • ■ A„(T), 7 (T)) satisfying 
(I3.42p . we obtain a unique hxed point of (I3.36p . Then 

k 

u{t, x) := U{t, a:) + ^ + 7 (t, x) 

2 = 1 

solves the equation fll.ll) on x [T, oo), satisfying 


iiv-kii 


L 7 >([T,oo);iflxL 2 ) < C6 


(3.43) 


with Lipschitz dependence on the data Aj(T) and {■y(T),'j(T)). Since the nonlinearity N only 
involves integer powers of 7 , we see that the integral terms in fl3.36p have smooth dependence 
on Ai, 7 . Hence we conclude that Xi{t),^{t,x) and the solution u{t,x) actually have smooth 
dependence on the data. 

By the estimates on A^ and 7 , we conclude in addition that 

n 

V = E^-(*)Pi +7((.7) e 7 C(|T,«d) X R3), 

2=1 

hence it (t) scatters to the same steady state as if (t) which is ( 0 , 0 ). 

We can now dehne 

T : T}{T) + (H,„((0,0)) nX,) L\ (3.44) 
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as follows: for any (70,71) e (^rid x 


and Aj G M such that 


e := A,(a, -hpi) + (70, 7 i) + ^(T) e f/(T) + (i?,„(( 0 , 0 )) nx,) , 

i=l 


set 

Ai(T) = Aj, for i = 1,... ,n and (7(T),7(T)) = (70,71)- 
Then with Ai(T) given by fl3.35p . we dehne 


^(0 := 


Xi{T)pi + 7o, Y HT)Pi 


7i 


, ^=1 


2=1 



If eo is chosen sufficiently small, then A, is uniquely determined by contraction mapping in 
the above. We dehne Ai as the graph of T and let XI be ^ {—T){M). We can then check 
that T, XI, XI verify the requirements of the theorem. Since (T) is a diffeomorphism, M. 
is a manifold. We remark that due to the presence of radiations the graph XI is in 
general not tangent to the center-stable subspace Xg. 

Step 3: unconditional uniqueness. Now suppose that we are given a solution u to equation 
(ll.ip . which satishes 

11 ”^ ~ '^llL°°([ 0 ,oo);hlxL 2 ) < Cl < eo. 

We need to show that lt{T) G XI. We denote 

n 

p{t,x) = u{t,x) -U{t,x) =Y + 7(^5 a;), 

2=1 

then if G L“([0, 00 ); x L^). By the fact that u, U are solutions to equation (II.ip and 
Strichartz estimates, we see that p G LfL^(J x M^) for any hnite interval I C [0, cxd) and 
admissible pair (g,r), we get for any T > T, 

ll^i(^)IU-([r,oo)) + ||7 (^)^)IIl 7> ([r,cx));hlxL2) <«i. (3.46) 

A.(i)ei"([r,r)). 

7 (i,i)ei 7 r(|r.r)xK=>). 


Notice the L°° bound on Aj implies that the stability condition fl3.35p must hold true, so we 
are again reduced to considering system fl3.36p . Now we wish to show that Xi{t) G 00 )) 

and l{t, x) G L‘fL’^{[T, 00 ) x M^). To do this, we follow similar arguments as in step 1. Dehne 
the norm 

n 

||(Al,--- ,Xn,j)\\x{[T,f)) ■ — ll^»(^)llL?(fT.rP + ll7llL?LS°('fr.f''lxR3~) 

2=1 
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By estimating fl3.36p similar to fl3.28p . fl3.29p . we get 


+ \\l\\LfL^([T,T)xR^) |-^i(^)l + II (7(^)) 7(^)) Ili/ixL^ 


+ 


2=1 2 = 1 

+ ll"^llLiL2([T,f)xIk“j ■ , 

Recall we have h = U — cf), ||h||i2^ooQ2-oo)xR3) < e = eo- Using the same estimate as in (I3.40p 
and (13.dip on the time interval [T,T), we obtain that 


II^IIl1L2([T,T)xR3) ~ '^o||hllL2ioo([r^f)xR3) W ^ 1 1 7 1 1 nL^^oo xRS) , 

k=2 


E 


and 


Hence 


< 


L^Lj ([T,oo)xII 


ll‘^ll7|ll^lll^=oL6([r,cx))xR3)ll''^ll7r7E([T’.«>)xR3) ~ ^1- 


i-\-j-\-k=b,k>l 


II (Ai, • • • , Xn, 7 ) llj\:([T,T')) ^ + ^o|| (Al, • • • , An, 7 ) |lx([T,T)) + ^ II (^1) ■ ■ ■ ) ^ni 7) II 


k 

x\T,fy 


k=2 


From this, by a continuity argument, we can conclude that 

II (Al,-- - , An,7)IU([T,oo)) < liminf ||(Ai, - - - , A^, 7 )||j,f([ 7 -j.)) < Uei < eo. 


T—>CXD 


and the contraction mapping theorem then implies 7t(T) G Ai. 

Step 4- summary. Let us sum up our construction as follows: consider any point (Uq, Ui) G 
At^, which generates a solution U(t,x) to equation (II.ip satisfying (I3.17p . For sufficiently 
large time T, we can construct a smooth graph Ai of co-dimension n in 

B.,( 7 (r)) e ky X 

such that solutions starting from Ai remain close to iJ (t) for alH > T and scatter to (0, 0). 
The graph can also be parametrized smoothly by 

Ai(T), - - - A„(T) G M, ^(T) G x 

in the following sense. For the parameters satisfying 

n 

+ 11(7(7), t(r))||*.,„<e„ 

2=1 

there exists a unique solution u to equation fll.ll) on t > T satisfying 

n 

u{T) = U{T) + E(T)p, + 7 (T, x), P^u{T) = ^(T) + P^dMT), 


2=1 


with the property that u(t) scatters to 0, and 

l7(i)-7(i)iL- ([r,oo),hixL2) 


<Ceo. 
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Moreover, any solution that satisfies ||^(t) —^(t)||^ixL 2 < ei with some ei < eo, for all times 
t>T necessarily starts on M. Using the solution fiow ^(t), we pull back our construction 
to time 0, ^(—T)A1, and the theorem is proved. □ 


Now we give a proof of Claim 13.21 Claim 13.21 will be proved as a consequence of the 
following lemma. 

Lemma 3.2. Let he a radial finite energy free radiation and (0, 0) be a steady state to 
equation jh. Recall that 

u = v/P^(-A - U + 504). 

Let 7 he the solution to 

dtt'j + oj'^'y = 0, in [T, oo) x 

^ (3-46) 

^(T) = P^(U^{T)). 

For any e > 0, if we take T = T(e, > 0 sufficiently large, then 

I L?L“([T,oo)xR3) < 


(3.47) 

Proof. For a given e > 0, fix 0 < 5 <C e to be determined below. We can take a radial smooth 




(3.48) 


compactly supported (in space) free radiation such that 

llU(o)-fo(o)||,j.,y„,, 

Let us assume that suppf/^(0) (1 Br{0) for some i? > 0. Hence by strong Huygens’ principle, 

for h 
that 


for large time T we have suppf/^(T) d Bt+r\Bt-r. Since f/^ is a free radiation, we see 


dttU^ - AU^ - VU^ + = -VU^ + in (0, cx)) x 


(3.49) 


By the decay property of V, 50"^ and the support property of U^, simple calculations show 
that 

tHSd ~ + 50'^f7^|lLiL2(['r^oo)xM3) = 0. 

Choose T sufficiently large, such that 

II — VU^ + 50'^U^||ii2,2([7^fx))xR3) < S. 

-A 


Note that if ^ — P^U^ solves 


(3.50) 


dttv + uj^v = —P-^ y—VU^ + , (t, x) G [T, cx)) x R'’^, 

with initial data ^(T) = P^ ^L^^(T) — P^U^fiP'^. By the bounds fl3.48p and fl3.50p . 

energy conservation for free radiation, and Strichartz estimates from Lemma 13.11 we can 
conclude that 

ll'^llLfLg°([r,oo)xR3) < C5. 


( 3 . 51 ) 
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Since f/^ is a finite energy free radiation, if we choose T sufficiently large, we have 

\W^\\l‘^L^{[T,qo)xR^) ^ C5. 


Combining bounds (I3.5ip and (13.521) . and hxing 6 small, the lemma is proved. 
Now the proof of Claim 13121 is easy. Note that due to the fact that 

||h(r) - (^,0) =0, 

we see that the initial data for 7 satishes 


(3.52) 

□ 


Jim || 7 (T)-P^f?yT)||, 


= 0 . 


T^oo. ' ..IHIXLW) 

Hence the claim follows from the above lemma and Strichartz estimates. 

4. Profile decomposition and channel of energy inequality 

In this section we recall some well-known properties of prohle decompositions hrst in¬ 
troduced in the context of wave equations by Bahouri, Gerard [3], and channel of energy 
inequalities discovered by Duyckaerts, Kenig and Merle [91CQ]. For both, we require the ver¬ 
sions adapted to the wave equation with a potential. We refer the reader to [16] for proofs. 
We hrst recall the following perturbation result. 


^ be an interval of time. Suppose u(t,x) G Ct{I, with 

a||^5/4^5/2^j^jg3^ < (3 < 00 and e{t,x), f{t,x) G LlLl{I x 


Lemma 4.1. Let 0 G / C 

ll'^llLfLiopxRS) < M < 00 , 
satisfy 

dttu — Ah -|- a{t, x)u + = e, (4.1) 

with initial data ^(0) = (ho,hi) G x Suppose for some sufficiently small positive 

e < Cq = eo(M, ft), 

ll|e| + |/|||lil|(/xr3) + IK'^OjMi) — (“0,hi)||/fixL2 < e. (4.2) 

Then there is a unique solution u E C{I, H^) with ||n||L5i^io(/xR3) < oo? satisfying the equation 

dttU — An -I- a(t, x)u + = f, (4.3) 

with initial data 5t(0) = 5t(0) = {uq,Ui). Moreover, we have the following estimate 


sup II^(t) 

t&i 




(^)llhixL 2 + ll"*^ ~ '^^IlnfLiopxRS) < C{M, ft)e. 


(4.4) 


Lemma 14.11 has the following implication concerning global existence and scattering for 
defocusing energy critical wave equation with potential, decaying both in space and time. 

Lemma 4.2. Let I be an interval of time and a G fl L\Lif{I x M^), and f G 

LjLl^I X Mf), with bounds ||a||^ 5 / 4 ^ 5/2 -|- ||a||i;^ii ;,3 < M and ||/||lil 2 < ft. Then there exists 

a unique solution u E C{I, H^) fl L^L)f{I x M^) to the equation 

dttU — An -|- a(t, x)u -f n^ = /, 


(4.5) 
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with initial data (mo,Mi) G x ril('^0)'*^i)|li7ixL2 ^ E). Moreover, we have 

ll'*^llLfLiO(/xR3) < C{E,M,(3). 

Thus if I = M, then there exist solutions uh, uf_ to free wave equation, such that 


lim ||M(t)-M+(t)||HixL 2 = 0, 

t^+oo 


lim ||u(t)= 0. 

c—>■—oo 


(4.6) 

(4.7) 

(4.8) 


The following lemma shows that for potentials a G (thus with space-time decay), 

large or small prohles are essentially not influenced by the potential. 


^IAt5I2. 


X 


and be a solution to the free wave equation 
Assume one of the following 


Lemma 4.3. Let a G L/ L 

m R X Take parameters {Xn,tn) with > 0, G 
conditions holds: 

1. tn = 0, lim (A„ 3 ^) = oo, 

2. lim ^ G {±cxo|. 

n^oo 

Let U he the nonlinear profile associated with U^, A„, tn- More precisely 

dttU -AU + U^ = 0mRx R^ 
with iJ (0) = (f/^(0), dtU^{A)) if tn = 0; or with 

.lim - ^WIIi7ixL2 = 0> (lim) 


£—^H “00 




(4.9) 


(4.10) 


if lim^^oo ^ = —oo (lim = oo respectively). Let Un be the solution to the Cauchy problem 


dttUn — Aun + a{t, x)un -|- = 0 in R X R^, 

™ttl?„(0)= (^C/‘(-^,|^),3;^9,C/'-(-^,^)). Then 

lim ( sup ||M^(t) - S(i)|liyixL2 + ll“n “ 

n^oo 

where f;), |;)). 


= 0 , 


(4.11) 


(4.12) 


The following prohle decomposition adapted for wave equation with potential plays an 
important role in our analysis. 

Lemma 4.4. Let a G L}J'^L}J'^ fl L(L^(R x R^). Suppose a radial sequence {uon,uin) G 
X is uniformly bounded and that we have the following linear profile decompositions 
(see Bahouri-Gerard m) 


j 

{Uon,Uin) = '^f(O) + 


j=2 '^jn 


L / tjn X 


_ _ _}_f)JjL(_'!2A Jf_ 

'^jn ‘^jn X ■ '^jTi '^jn 

jn 


)) -|- u)’j„(0), (4.13) 
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with the following properties: 


Uj, wjn are radial and solve the free wave equation for each j, J; 


either tjn G M, Xjn > 0, lim G {± 00 } or tjn = 0, lim ( Aj„ + 


jn 


A 


•jn 


= OO] 


r ■ / ■/ 1- f I ^j'n , \l'jn tj'i , , 

for J / j', hm + -f -h —— I = 00 ; 

n^oo \Ajin Ajn Xjn 

write Wjn{t,x) = then 0, and wjn ^0, as n ^ oo] 

jn 

lim lim sup ||wj„|L5iio(RxK3) = 0. 

J—^CO Yl - ^QQ 

Let Ui satisfy 

dttUi - AUi + a{t, x)Ui + Uf = 0, inR x R^, (4.14) 

with f^i(O) = 17f(0). Let Uj be the nonlinear profile associated to U^, Xjn, tjn as defined in 
Lemma \4 . 3\ for j > 2. Let Un he the solution to 

dttUn — AUn + a{t, x)Un + = 0, iuR X R^, (4-15) 

with t^n(O) = {uQn,Uin)- Then we have the following decomposition: 

J 

= lfi{t) + n{t) + ^jn(t) + Jn{t), 


(4.16) 


1=2 


with 


lim limsup ( sup ||^Jn(i)||i/ixL 2 + lkjn||LfLiO(RxR 3 ) j = 0, (4.17) 

J-^oo n^oo Vi 6 K / 


where Tfjn{t,x) = ( ■:^Uj(*-^, ^),— 


^L/2^j\ ) \. )l , 3 / 20 'tOjq , ’A 

Ad 'Ajn Ad 'Ajn 'Ajn 

jn jn 


Moreover, denoting Uin = Ui, 
for pn > (Tn > t) and 9n eM we have the following orthogonality property for 1 < j ^ j' 

' (4.18) 


lim / ^Ujn^Ujfn “t“ d)tUjnhdtUj^nidn, dx 0, 

n^oo Jcrn<\x\<pn 

lim / S/UjnS/Wjn + dtUjndtWjn{0n, x) dx = 0. 

Ja-n<\x\<p„ 


(4.19) 


We also need the following channel of energy inequality from [TB] , which was proved with 
similar arguments as in Duyckaerts, Kenig and Merle HD]. 

Theorem 4.1. Suppose radial finite energy {uo,Ui) ^ (0,0) for any steady state solution 
(0,0) of eguation Let n G C(R, 77^) fl LjL^°((—T, T) x R^) /or any T G (0, cxo) he the 

unigue solution to eguation M-l\) with t^(0) = {uo,Ui). Then there exists R> 0 and 5 > 0 
such that 

[I Vm|^ + {dtu)'^]{t, x) dx > 6 > 0, (4.20) 


' p|>r?+|t| 


for allt>0 or all t < 0. 
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This theorem tells us that if a radial solution to equation fll.ip is not a steady state, then 
it must emit energy to spatial inhnity. For applications below we also need the following 
quantitative version of Theorem 14.11 


Theorem 4.2. Take 1/ G hi C T. Suppose that radial finite energy initial data {uo,ui) ^ 
(0,0) for any steady state solution of equation U.l\) . with ||(mo, wi)||^ 1 x£, 2 (k 3 ) < M < oo. Let 

u G C{R, H^) n 4^L^°((-T, T) X M^) 

for any T G (0, oo) be the unique solution to equation M.l\) with ^(0) = {uo,Ui). Denote by 
S the set of radial steady states of equation M.l\) and define 

6 := inf {||(mo,Mi) - (0 ,O)||^1 xl2 : (0,0) G E} > 0. (4.21) 

Then there exists c = c{6, M,V) > 0 such that 

f [\Vu\‘^ + {dtu)‘^]{t,x) dx > c, (4.22) 

J\x\>\t\ 

for all t >0 or all t < 0. 


Proof. Suppose the theorem fails. Then there exists a sequence of solutions Un to equation 
([HI]) with initial data {uon,Uin) e x satisfying ||(Mon,Min)|lij-ixL2(M3) < M, (14.211) 

with a uniform 5 > 0, and 

inf / \Vun\'^ + {dtUn)'^it,x) dx + inf [ \Vun\‘^ + {dtUn)'^it,x) dx <-. (4.23) 

By passing to a subsequence we can assume that (won^Wi^) admits the following prohle 
decomposition 


ip^Ori) ^In) 



P) 1 + ./»(•)). (4-24) 


X 


Xin X 


■jn ''jn 


X 


3/2 

jn 


Xjn X 


jn 


with the following properties: 


Uj, wjn are radial and solve the free wave equation for each j, J; 
either tjn G M, Xjn > 0, lim G {±C)o} or tjn = 0, lim [ Xjn + t— ) = oo; 


A, 


n—>oo \ 

A,- 


jn 


A 


•jn 




forjV/, lim ( Till + Tin + 1''^- 


Xj'n Xjn 


X 


'jn 


= oo\ 


write Wjn(t,x) = -^)y then 0, and wjn ^0, as n co; 

Xj!^ ^jn ^jn 

lim lim sup ||u;j„||i5iio(RxK3) = 0. 

J—^CO - ^QQ 




< OO, 


Note that 
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thus by finite speed of propagation we can apply Lemma 14.41 in the exterior light cone 
{{x,t) : |a:| > |t|}, and obtain 


j 

= + ^T}jn{t)+ + for |a:| > |t| (4.25) 


i=2 


with 


lim limsup ( sup |l^Jn(t)||i/ixL 2 + |kjn||LfLiO(RxR 3 ) ) = 0, (4.26) 

where iJjnit, x) = i , ^), ^) ], and with Uj and Ui given as in 

_ \^jn _ / 

Lemma US] and Lemma 031 Moreover, denoting Uin = Ifi, for p„ > cr„ > 0 and G M we 
have the following orthogonality property for 1 < j ^ j' 


lim / ^ LIjfi^ Uj'fi dtUj'fidtU>]ri^ x) dx 0, 

J an<\x\<pn 


lim 


'a„<\x\<pn 


^Ujfi^wjfi dtUjffdtwjfi{0fi^ x^ dx 0 . 


(4.27) 

(4.28) 


If for some 2 < j < J, ^ 0, by results in [16] (see remark at the end of proof of Lemma 4.6 


m 


there exists some fixed e > 0 such that 


hl>l4 


iVt/jnl^ + {dtUjnf{t,x)dx > e > 0 , 


(4.29) 


for alH > 0 or all t < 0. Thus by the orthogonality property of profiles fl4.27p . fl4.28p we have 


inf / 

> - inf [ 

2i>0^|>|, 

> - > 0 
“ 2 


IVUnl"^ + {dtUnY{t,x)dx + ini / \'VUn\‘^ + {dtUnY{t, x)dx 


^^^J\x\>\t\ 


IVfojnp + {dtUjnY{t, x)dx + - inf / \VUjn? + {dtUjnf{t, x)dx 


2 J\x\>\t\ 


for all n snfficiently large. A contradiction with fl4.23p . Thus we must have Ly = 0 for 
2 <J< J. The profile decompositions of (unOj^ni) then simplify to this form: 

(M 0 n,Mln) = '^f(O) + Mt„(0). (4.30) 


By the decomposition fl4.25p . Theorem 14 .1 1 and orthogonality property of profiles, using the 
same arguments as above, we conclude that ^f(O) must be a steady state. Thus from the 
bound fl4.2ip we have 

ll^n(0)|| hixL2(R3) — 


( 4 . 31 ) 
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By the channel of energy estimates for the linear wave equation, by the decomposition fl4.25p . 
and the orthogonality property, we obtain 



for all sufficiently large n. We thus again arrive at a contradiction with fl4.23l) . The theorem 


is proved. 


□ 


5. Center Stable manifold of unstable excited states 


In this section we show that any unstable excited state can only attract a hnite co- 
dimensional manifold of solutions and hnish the proof of our main Theorem 11.31 We hrst 
prove Theorem 15.11 which provides the key estimate on the energy of the radiation term. 

Before going into the technical details let us briefly outline the main ideas underlying 
the proof. By results in Section [3l we know that if a solution U{t) scatters to an unstable 
excited state (0, 0), then there exists a local hnite co-dimensional center-stable manifold 

around ^(0), on which solutions scatter to the same excited state. We would like to show 
the following in a small neighborhood of if (0): if the data {uq, ui) does not he on this local 
center-stable manifold, then the solution u (t) with initial data {uo,ui) will scatter to a 
steady state with strictly less energy, thus not to (0, 0). Then it is clear that the set of initial 
data in x for which solutions scatter to an unstable excited state is a global hnite 

co-dimensional manifold. 

Note that the local center-stable manifold theorem 13.21 guarantees that the solution u{t) 
will exit a small ball centered at ^(0). However after it exits the small ball, we will lose 
control on the solution based on perturbative analysis alone. Thus we need some global 
information about the future development of the solution u{t). The key global information 
here is the channel of energy inequality. Roughly speaking we show by a channel of energy 
argument that u{t) will emit at least some hxed amount of energy more than solution starting 
from the center-stable manifold to spatial inhnity, thus leaving u{t) with less energy in a 
bounded region as t —)■ cxo than is required for it to scatter to (0,0). This forces u{t) to 
scatter to a diherent steady state from (0, 0). The precise result is as follows. 

Theorem 5.1. Let V E . Suppose that radial finite energy solution 'if {t) to eguation 

(toy scatters to an unstable excited state (0,0). Let M. he the local center-stable manifold 

around ^(0) and let eo, ei be as defined in Theorem. IS. SI Then there exist e with 0 < e < 
ei < Co and 5(ei) S> e, such that for any solution u with radial finite energy initial data 
{uq, ui) ^ Ai with 
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f 

\\Vut , 

J\x\>t—L 

2 2 


we can find L > 0 such that 

{t,x) dx > S(lf (t)) — S{{(j),0)) + 6, for t > L. (5.1) 

Suppose St(t) scatters to (0i,O) G S (the set of steady states). Then 

£((01,0)) <£((0,0)). (5.2) 

Proof. By the local center-stable manifold theorem of Section [3l the locally dehned hnite co- 
dimensional manifold M satishes the property that any solntion to eqnation fll.ll) with initial 
data on M. scatters to (0, 0). Moreover, if a solntion it (t) with initial data (mq, ui) G (jJ o) 
satishes 

\\lt{t) — it)\\mxL2 < <^1 foi’ all f > 0, (5.3) 

then (mo,mi) G Ai. By shrinking ei if necessary we can assnme that the distance from any 
other steady state to (0,0) is greater than 3ei. Take e < ei snfhciently small to be chosen 
below. Since solntion U (t) scatters to (0, 0) as f —)■ cxo, denoting by the scattered linear 
wave, we have the property that 

lim 10(t) - Jf^it) - (0,O)||j 


£—>-cx) 


\mxL^ 


= 0 . 


(5.4) 


By 05.41) . the fact that 0 G and G LfL,),°([0, cxo) x M^), for any small (5i > 0, we 

can choose L > 0 and Ti> L snfhciently large snch that for t >Ti we have: 

• (Closeness of if to + (0, 0) and choice of the bonnded region) 

\0{t) — lj^{t) — ( 0 , 0 )||^lx ^2 + ||0|lii-l(|3;|>L) — (5-5) 

• (Most energy of the free radiation is exterior) 

[ \17x,tU^t{t,x)dx> [ \17t,xU^t{t,x)dx - 5i; (5.6) 

J\x\>t-Ti+L 0M3 

• (Control on the Strichartz norm of the radiation) Let 

D := {{x,t) : \x\ < Ti + L — t, 0 < t < Ti}. 

Then we have 

l|l^^llLfLiO((0,oo)xR3\D) < (5.7) 

We remark that 05.61) ensnres that can essentially be taken as zero for onr pnrposes inside 
the region |a;| < f — Ti -|- L for t > Ti, which will be important to keep in mind later, in 
order to distingnish the second piece of radiation. 

By the continnons dependence of the solntion to eqnation 01.ip with respect to the initial 
data in x L^(M^) and by hnite speed of propagation, if we take e snfhciently small and 
radial initial data {uo,Ui) G x Lf\Ai with ||(mo,Wi) — < e, then 


(5.8) 
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can be made sufficiently small. Hence, noting that ||H|L5/4r5/2,| i,,. is finite, we can apply 

Lemma [4. II to conclude that 

-Tf {t)\\HixLmx\>t-T,) < ^ 1 ) for all t > Ti. (5.9) 

fl5.9p means that we can effectively identify it with iJ in the exterior region 

|x| > t — Ti, t > Ti. 

Hence by fl5.5l) . we see 

\\lt{t) — ij {t)\\H'^xL'^(\x\>t-T-i_+L) — (5.10) 

that is, we can also identify it with in the exterior region |a:| > t — Ti + L, t > Ti. 

At this point to avoid any possibility of confusion due to the many parameters, we remark 
that and e can be made as small as we wish, with Ti, L depending on and iJ only, e is 
a small free parameter below some threshold determined by ^i. The key point for us is that 
ei > 0 is fixed no matter how small e is chosen. 

Since {uq^Ui) ^ A4, there exists T 2 > 0 such that 

||V(r2)-h( ^2)||i/lxL2(R3) — ei. (5.11) 

Note that the choice of Ti and L does not depend on e, thus by the continuous dependence 
of solution on initial data in H 1 X L2, if we choose e sufficiently small, we can assume 
T 2 > L + Ti + 1. Let us consider the data lt(T 2 ) in more detail. By estimates fl5.5p and 
fl5.1ip we can write 

lt{T2) = {<!),0) + l}\T2) + lt, 

where it G x satisfies 


(5.12) 

(5.13) 


2 ei > ei + (5i > ||^||i/ixL2(iR3) > ei - > ei/2, 

if 5i is chosen smaller than y. Consider the solution u{t) to equation fll.ip with 

'u{T2) = ( 0 , 0 ) + it. 

Then by the quantitative channel of energy inequality from Theorem 14.21 we infer that 

/ \17t,xutit, x)dx > c(ei) > 0, for all t > T 2 or all t < T 2 . (5-14) 

J\x\>\t—T 2 \ 

By bound fl5.7l) and Lemma [4.11 we obtain that for |a;| > |f — T 2 I 

it (f, x) = x) + ^{t, x) + 1^ (f, x), (5.15) 


where the remainder term 1^ satisfies 


SUp||^(f)||^lxi2(R3) 
iGR 


< CSi. 


(5.16) 
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We claim that the channel of energy inequality fl5.14p holds for all t > T 2 . Otherwise, 
inequality (I5.14p holds for all t < T 2 . By (I5.15p and (15.161) . setting f = 0 and noting that 
\W^\\h^xl‘^(\x\>T 2 ) made smaller than S if T 2 is large enough, we see that 

11 (^0, Wl) 11^1 > c(ei) - > ^c(ei), (5.17) 

if is chosen sufficiently small, a contradiction with hniteness of t/(0) in x and 

IKwo^wi) — ^(0)||^ixL2 < e for T 2 large. Thus we have the following channel of energy 
inequality 

/ dx > c(ei) > 0, for all t > T 2 . (5.18) 

J\x\>t—T2 

The estimate 05.101) . decomposition 05.151) and estimate on the remainder term 05.16P imply 
for f > T 2 

I Vi_a;'Up(f, x) dx < C6i, 


\x\>t—Ti-\-L 


and consequently 


Hence 


f t — T\-\-L'^\x\'^t —T2 


^ C( 6 i) - 


(5,19) 


(5.20) 




'H^xL'^{\x\>t-T2) 

^ ll'^^llitixr2(K3) + c(o) - C5i, 


if is chosen sufficiently small. Choose (5i sufficiently small depending on ei, 05.ip is then 
proved with some 5 = (5(ei) > 0. 

Now let be the free radiation term for the solution lt{t) as f ^ 00 . Then 

ll^(i) - - (<(>1,0)11 ifixL2(R3) — 0. (5.21) 

By the first part of Theorem 15.11 


II “ WmxLHm - 


> lim ||5t|| 


> 

irixL2(|x|>i-r2) — 


Il 7 


L ||2 

hixL2(R3) 


+ c{ei). 


(5.22) 


Note that 


lit? 


L ||2 

H^xL'^ 

I WmxL^ 


f(I?)-£(.)., 0)i 

= S{lt) - £?(01,O). 

If ||'i?(0) — ^(0)||^ixL2(R3) is chosen small, by 05.221) we conclude 
The theorem is proved. 


(5.23) 

□ 





















35 


Proof of Theorem \1.3[ We only consider the case in which (0,0) is unstable, the case of (0, 0) 
being stable can be handled using standard perturbation arguments. Since in some small 
neighborhood of any point iJ ( 0 ) on coincides with the local center-stable manifold 

Ai of codimension n which we constructed in Section [3] by Theorem 15.11 is thus a global 
manifold of co-dimension n. The path-connectedness follows from the following theorem. □ 

Theorem 5.2. For any unstable excited state (0,0), the corresponding center-stable manifold 
is path connected. 

Proof. Given data (mq, wi), (uo, wi) £ we denote the corresponding solutions by u,u. 
Write h = u — (f),i = it — (j). Repeat step 0 and step 1 in the proof of Theorem 13.21 Then 
given any e 1, we can hnd T = T(e, u, it), such that 

\\h\\L^L'^nL^L%([T,oo)xR3), ||^||L2L“nLf=L|([r,cx))xR3) < (5-24) 

Now we seek a function w{6, t, x) of the form 
w{6, t, x) = (1 — 6)u -i 6u -if] 

n 

= 0 -F (1 - 6)h + 6i + ^ Xi{6, t)pi + 7 ( 0 , t, x) (5.25) 

i=\ 

such that for all 6 G [ 0 , 1 ], 7 ( 6 ^, t,x) P pi,i = 1,... ,n and w{6, t, x) is a solution to equation 
fll.ip that scatters to (f. 

For 9 G [0,1] hxed, the equation satished by 7 = + 7(^5 i®- 

ptt - Ap -V{x)p + -F N{6, h, i, 0, p) = 0, 

where 

N{e, h,i, 0, 7 ) = (0 + (1 - e)h + ei + p)^ - (i - 0)(0 + h)^ - 0(0 + - 50 ^ 7 . 

Now we can repeat the stability condition fl3.35l) and obtain the reduced system of the form 

In N{d, h, I, 0, p), the terms independent of p are of the form 

{(f + {1 - e)h + ee)^ - {1 - e){(p + h)^ - e{(p + i)^ = c{9,i,j,k)fhH’^. 

i+j+k=5,i<3 

Notice that there are no terms 0^ or cf^h, 

Also, the linear term of p in N{6, h, I, 0, p) is 

5(0 + (1 — 9)h P 9t)^p — 5(j)‘^p 

hence all linear terms involve a factor of h or i. 

Now we can repeat estimates 03.381) fl3.39p . then O3.40p for the linear term in p, 03.411) for 
higher order terms in p. We also have the following estimate on terms independent of p 

Y, c(e,i,j,k)4,‘hir 

i-\-j-\-k=b,i<3 



LiL2([T,oo)xR3) 
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For example, one checks that 

I|0^^^IIl1l2([T,oo)xR 3) < ll^lliell^llisj^go ^ 

using (I5.24p . To sum up, using the X norm dehned in (13.371) . we conclude that 

IKA,.-- - ,A„,7)IU(|T,co))<A't= + A'(X;|Ai(9.r)| + ||(7(9,r),7(9.r))||,j.,p 

\i=l 

5 

+ -^'i^e||(Ai, • • • , A„,7)||x([r,cx))) + K 

k=2 

where K is some absolute constant. 

Moreover, in a similar fashion one sees that the difference of two solutions satishes a similar 
estimate in which the hrst two terms disappear. Using the contraction mapping principle, 
we conclude that for sufficiently small data 

n 

\\{0,T)\ + ||(7(«.T),mr))|ls.,„ < i 

i=l 

there is a solution w as in fl5.25l) which solves fll.ip . We can also check that w scatters to 0. 

In particular, let us take Xi{9,T) = ^59(1 — 9) and j{9,T,x) = 0. We claim that the 
corresponding solution w{9,t,x) satishes the following relation 

w{0,t, x) = u{t, x), x) = u{t, x), for alH G M. (5.26) 

In fact, notice that Ai(0,T) = 0,7(0,T, x) = 0 implies Ai(0,f) = 0,7(0,t,a:) = 0 for f > T, 
which further implies t(;(0, t, x) = u{t, x),t >T. Similarly we have w{l, t, x) = u{t, x),t >T. 
Then fl5.26l) follows from the uniqueness of solutions to equation fll.ip . 

Hence {w{9, 0,x),9 G [0,1]} is a path in connecting the two data {uq, ui), {uq, hi). □ 

6. Appendix A: Some elliptic estimates. 

We begin with the following lemma. 

Lemma 6.1. Denote by Bi the ball of radius 1 in M^. Let V G and A > 0. Suppose 

u G H^r\L^{Bi) is a weak solution to 

— Am — Um + Am® = 0, in Hi, (6.1) 

in the sense of distributions. Then u G L°°(Hi) and 

||m||l°°(Bi) < C{\\V\\L^(^Bi))\\u\\L^iBi)- (6.2) 

Remark 6.1. The assumption on the regularity of V can be signihcantly relaxed to U G 
with g > |. For the sake of simplicity, we shall not prove the most general version. The 
proof is based on DeGiorgi-Nash iteration arguments. 
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Proof. Assume u is not identically zero. By multiplying a positive constant to u, we can 
assume in addition that = 1. u still satisfies equation (16.Ih with possibly a different 

A in the equation. Fixing an M > 1 to be determined later, we need to show ||n||ioo( 5 j) < M. 

For each integer A; > 1, we define 

1 1 

c, = {l-^)M. 

For each /c > 2, we fix a smooth nonnegative cutoff function such that 

supp? 7 fc (E , with 0 < < 1 and iVr^^l < C2^. By a sign change, it suffices to prove 

2 

that u{x) < M ioT X ^ Bi. Multiplying equation fib.ip with {u — Ck)+Vk and integrating, we 
get that 


V(m - Cfc)+V ((u - Ck)+r]k) dx- V{u- Ck)+ ■ {u - Ck)+r]kdx 




Bi 

5 / 


CkV {u - Ck)+r]k dx + X u^{u - Ck)+r]k dx = 0, 


>Bi 


iBi 


which implies that 


Thus, 


|V(m - Cfc)+p + (n - Ck)+dx 


'B, 


'^k 


< C{\\V\\loo(^Bi) + 1 ) 


'B 


{u — Cfc)^ dx + 


^k+'^k-i 

2 - 


IB 


{u — Ck)+ dx + C4f / [u — Ck)\ dx. 


'^k+'^k-i 

- 2 - 


IB 


2 - 


’ Br 


{u - Ck)% dx 

<cm\L-iB,) + i)^^ 

-\-\\y\\L°°Ck 


{u — Ck-i)\ dx + 


' B. 


'^k-1 


/ {u-Ck)\dx 

2 

{x e Br^+r^_i : u{x) > Cfcl 

\ d Br^+rf._^ j 


2 


Note that by Chebyshev inequality we have 


{x e B r^+r^_^ : u{x) > Ck] 


< 


(Cfc Ck—\) 


' B.. 


{u - Ck-i)\dx. 


^k-1 
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Hence, we get that 


Note that 


{u - Ck)% dx 


'B, 


'^k 


C/c j 


(m - Ck-i)\dx + 




^k-1 


+ll^lli 


Cfc ^k—l JBr 


<C'(II^IU^(b,) + 1)4^ 


{u - Ck-i)+dx 
{u - Ck-i)1dx. 


' Br 


{u - Ck)% dx> {u- Ck+i)\ ■ (cfc+i - Ckf dx. 


’ Br 


'B., 


^k+i 


Summarizing the above inequalities, we obtain 


{u-Ck+iY+dx] < C(||H||lcx,(b^) + 1)2“^^ / {u-Ck-i)Xdx, 


fB. 


'^k+i 


^ Br 


Denote 



then equation (16.3p can be written as 



6,+i < C^16>^eU < (C'il6)"6ti, 


for all k >2. 

(6.3) 

(6.4) 


(6.5) 


where we have suppressed the dependence of Ci > 1 on H. Recall that a routine energy 

inequality implies ||Vn||L 2 (B 3 ) < C||M||i;, 2 (B 3 ) = C. Hence, if we choose M sufficiently large, 

% 

by the Holder and Chebyshev inequalities, we can assume e 2 , es satisfy 


^2; ^3 ^ C —r < 


Mi 


(48Ci 


( 6 . 6 ) 


Then from the iterative inequality fl6.5p we can prove by induction that for all 

k > 2. Hence we have 


lim 

k—^oo 



Ck)+ dx 


0 . 


(6.7) 


Note that | and Ck ^ M as k ^ oo, thus we conclude u{x) < M ior x E Bi. The 

lemma is proved. □ 


Now we are ready to prove the main result on the regularity and decay of steady states to 
equation fll.ip . 
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Theorem 6.1. Let V E Y. Suppose u G fl is a steady state solution to 

equation t.e., u solves 


-Au-Vu + u^ = 0 in 

in the sense of distributions. Then u G (M^) for any 1 < p < oo, and 


\u{x)\ < 


C 


(1 + kl)’ 


X G 


( 6 . 8 ) 


(6,9) 


Proof. By Lemma f6.1l we see u G Then the estimates follow from standard 

elliptic regnlarity theory. Let ns now tnrn to the proof of the decay estimate fl6.9p . For any 
2R = l^ol > 2, set v{x) := R^u{Rx + Xq). Then v solves 

-Am - Vr^xqU + m'’ = 0, in i?i, 

where Vr,xo = R^V{Rx + xo). Since V eY with snp(l + |x|)^|l/(x)| < cxo for some (d > 2, 
we see that ||Vr,xo|U°°(Bi) < ||l^||r- Thus by Lemma [6Tl 

|' y ( 0 )| < \\v\\l2(Bi) < C\\v\\l6^Bi) < C\\u\\l6^{i^i>b/2}). 

Hence by rescaling and the fact that Xq is arbitrary, we see that 

1 


|m(x)| = o 


(1 + |x |)2 


as X —)■ oo. 


Now choose i? > 1 sufficiently large, and set m(x) = R 2 u{Rx), then vIqb^ = R^uIobj^ is 
small, and Vr := R?V{Rx) is small in the exterior of i?i, in the sense that 

sup 1(1 + \x\YVr{x) \ < ei{ —)■ 0, as R ^ oo. 

\x\>l 

Then for sufficiently large R, by a standard perturbation argument one constructs a solution 
V to the equation 

— Am — 1/rm + = 0 , in (6.10) 

with the boundary condition v\dB-i = with the estimates 

cr 


|m(x)| < 


(1 + kl)’ 


for |x| > 1. 


( 6 , 11 ) 


Since both v and v are small solutions in to equation fl6.10p with the same boundary 
condition, we conclude v = v. Hence v also satishes (16.111) . Scaling back to m, we see that u 
satishes 

C 

\u{x)\ < - ----, for X G (6-12) 


(1 + kl)’ 


The theorem is proved. 


□ 
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7. Appendix B: Endpoint Radial Strichartz Estimate 

In this appendix, we give a proof of Theorem 13.11 in the radial setting. We only need 
to prove the endpoint case, since other cases are known from [13]. In the case of the ho¬ 
mogeneous equation, Klainerman and Machedon [20] first observed that the endpoint case 
(g, r) = (2, cx)) holds true for data of the form (n(0),ni(0)) = (0,g) for radial function g. 
Here we extend the result to the inhomogeneous equation with general radial data. 


Theorem 7.1. Let v be a finite energy solution to the 3d wave equation 

{dtt - A)v = F 

with initial data (n(fo), 5 r(A)) = {f,9) G 77^ x and assume that f,g,F are radially 

symmetric. Then we have the estimate 


ll«llLfL-(Ex..)<C(«,r)(||(/,j)|| RixL2 + \\F\\ LiLi(RxR3)) , 

Proof. The solution takes the form 


v{t, x) = cos(| V|f)/ -h sin(| V|f)g -h 


sin(|V|(t-s)) 


F{s)ds. 


(7.1) 


Case 1: / = 0, F = 0, this was proved in [20] using the explicit formula 


v{t, r) = 


ft+r 


'\t-r\ 


9{p)pdp, 


t > 0. 


Denote M{-) to be the Id maximal function defined as 

M{f){x) = sup^ [ 

xei b I Jj 


I/I 


for any interval J C M. Extending g to be an even function defined on M, we then have 

^ rt+r 


sup \v\ < C sup 

r>0 r>0 


\i-A 


\9{p)p\dp 


< CM[g{p)p\{t). 


This is obvious for r < t, while for r > t, notice that [t — r,t + r] is a larger interval than 
[r — t,t + r], hence the term in the middle is still bounded by maximal function. 

Using the Hardy-Littlewood maximal inequality [23], we get that 


sup |u|||l 2 (m+) < C\\M[g{p)p]{t)\\L 2 (^+) < C\\g{p)p\\L 2 (R) = CyWiHRS). 


r>0 


By time reversibility, we have the same estimate when integrating on t G (—oo, 0]. 
We rewrite the estimate in the following form 


■sin(|V|t)g 


< 


L^L“(RxI 


I|5||l2(E3)- 


(7.2) 
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Case 2: f = g = 0. Using fl7.2p and Minkowski inequality, we immediately get 


sin(|V|(t-s)) 


F{s)ds 


< / 

L?Lg°(RxK3) >''0 


sin(|V|(t-s)) 
xm (s)-- F{s) 


ds 


(Rxl 


||-^('S)||l2(R3)C^S — ||-^||l1L2(RxR3)- 


Case 3: = 0, F = 0. For simplicity, we prove our estimates for Schwartz function /, and 

the estimates for general / G H^ad follow by approximation. 

In the radial case, we have the following explicit expression for v: 


i>r+t 


v = dt- f{p)pdt =-[f{r+ t){r+ t) - f{\r-t\){t-r)], t > 0. (7.3) 


As a first step we bound 


sup \v{t,r)\ < sup - 

0<r<t 0<r<t f' 


ft+r 


ipfip)ydp 


' t—r 


M[{pf{p)y]{t). 


Applying Hardy-Littlewood maximal inequality and Hardy’s inequality [23], we get 

^(P 

/ 


sup |u|||l2(r+) <||/(p)+ p/'(p)||l2(ir+) 

0<r<t 


< 


\x\ 


1,2 (]R3) 


+ l|V/||t>(..) < ||V/|U=,K.). (7.4) 


Next if r > t we claim that 


sup|/(p)|||i2(K+) < / |/'(p)|Vdp. 

P>t ' Jo 


(7,5) 


Dualizing 07.51) we see that it is equivalent to 


f'{w) dw h(t, p) dp dt 


p>t>o J p 


wf{w) 


h{t, p) dpdt 


w>p>t>0 


dw 


< 

r\j 


II/'(«;)«; ||l2^(r+) II h||i2M(R+xR). 


(7.6) 


So we need to prove that 

1 


w 


h(t, p) dpdt 


w>p>t>0 


L2(R+) 


,7) I|U||l2j;^1(M+xR)) 


by change of variables t = wr, we have 

\h{wT, p)| dp dr 


W>p>WT>0 


< 


L2(R+) 


1 /‘OO 



0 JO 


< 

r\j 


< 


\h{wT, p)| dp dr 
\h{t,p)\dp 


Ll 


dr 

l2(e+) 


r\j 











































as we claimed. Thus we proved fl7.5p whence 


Next we will prove that 


sup -|/(r + t)\{r + t) 

r>t >0 r 


sup -\f{r-t)\{r-t) 

r>t>0 f 


< 


L2(K+) 


llv/IU- 


< 


L2(R+) 


IIV/IU 2 . 


By the change of variables r = t{p + 1), this is equivalent to 


sup 


P 


p >0 1 + P 

By duality, this is equivalent to 


I/Ml 


< 


r*00 \ ^ 

l/'(p)iv^p • 


poo poo 


p 


or 


'0 Jo 


wf{w) 


1 + p 


f'{w)dw h(t,p) dpdt 


tp 


< 


||/'(-u;)-u;||l 2 (R+)||h||i 2 M(R+xR), 


P 


w>tp >0 


1 +p 


So we only need to show that 


P 


Iw>tp>o ^ T p 
By the change of variables t = tw 


h(t, p) dpdt 


h(t, p) dpdt 


dw 


||/'(w^)w^|U2(R+)||h||i2M(R+ 


< 


Ll 


LlL+R+xl 


P 


w>tp ^ T p 

P 

l>Tp 1 + P 


h(t, p) dpdt 


h{TW, p) dp dr 


L2(R+) 

Ll(++) 


< 


< 

rs_/ 


< 


'1 Jo 

“1 roo 



0 JO 


|h(rtc, p)| dr dp 

\h{TW, p)| dpdr 
dr 


+ 





p\h{TW, p)| dr dp 


0 Jo 


Ll 



POO ^min(i,l) 

+ 

/ / p h(rt<;,p) dpdr 

L2(R+) 

Jo Jo 


Ll 




<11/7,11 r2 


L 2 l 1 (R+xK) + 


_ + 

00 2 


^00 poo 


'0 Jo 


l + r 


\h{TW, p)| dpdr 


Ll{R+) 


'0 


l + r 


L 2dr||h||i2^i(R+, 


< 


rsj 
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Hence combining fl7.3p and 07.71) . 07.81) . we deduce 

l|sup|n||li2(]R+) < ||V/|1 l2(jr3). 

r>t 

Together with 07.41) we have proved 

lkllL2L“(R+xR3) ^ ||V/||l2(r3). 

By time reversibility, we get 

lkllL2Loo(RxR3) < ||V/||l2(r3) 

and we are done. □ 

Acknowledgement: The proof of the uniqueness of ground states up to a change of sign was 
communicated to us by Tianling Jin. 
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